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FINAL REPORT ON 


MEASUREMENT OP THE VISCOSITY BY THE OSCILLATING CRYSTAL METHOD 
OF FOUR FLUIDS AS FUNCTIONS OF PRESSURE, 

SHEAR RATE AND TEMPERATURE 


by R. G. Rein, T. T. Charng, C. M. Sliepcevich, and W. J. Ewbank 
UNIVERSITY OP OKLAHOMA RESEARCH" INSTITUTE 

INTRODUCTION 

This report summarizes the work performed under Contract 

No. NAS 3-13487. The viscosities of four fluids supplied by 

NASA were measured over a range of pressures and temperatures 

using an oscillating quartz crystal method as described by 

Rein, et al . (13). The four fluids supplied by NASA were; 

1. Mobil XRM-109P, a synthetic paraffinic hydrocarbon. 

2. DuPont Krytox 143-AB, a perf luorinated polymer. 

3. Humble FN-3158, a super-refined naphthenic mineral oil. 

4. Humble FN-‘3158 fluid plus ten percent volume of Kendall 
0839 super-refined paraffinic heavy resin. 

The following test conditions were specified for each fluid: 

1. At a constant temperature of 100 ®F and a constant pressure 
of 10,000 psig using crystal frequencies of about 20 kHz 
and 60 kHz. 

2. At a constant crystal frequency of about 20 kHz and at 
pressures of 0, 10,000, 20,000, 40,000, 60,000, 80,000, 
100,000 and 115,000 psig at constant temperatures of 100, 

210 and 300®F and 250®C (482®F) . 
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The current limitations on equipment restricted the 
measurements of viscosity on the four test fluids to a pressure 
of 40,000 psig and to a temperature of 300 ®F. However, it 
should be noted that density measurements were made up to 482®F 
and pressures to 115,000 psig with the calibration fluid di- 
(2-ethylhexyl) sebacate, which demonstrates that the present 
equipment is operable over the full temperature and pressure 
range originally sought but not for very viscous (above 10 
poise) fluids without equipment modifications. 

The following sections in this report include a descrip- 
tion of the experimental techniques, the results in graphical 
form and tabulated experimental data in Appendices A and B, 
conclusions, experimental limitations and recommendations. 



EXPERIMENTAL TECHNIQUES 


The measurement of viscosity of fluids as a function of 
pressure has been reported by a ntomber of investigators, but 
the most comprehensive and probably best known early studies 
are those of Bridgman (5) on pure compounds and the ASME 
Pressure-Viscosity Report (2) on lubricating oils. In both 
instances a falling-weight technique was used in which the 
shear rates were far too low to detect non-Newtonian behavior. 
Most of the published experimental work to date on viscoelastic 
behavior of fluids has been at atmospheric pressure. 

The first significant studies on non-Newtonian behavior 
of fluids (polymer solutions) at high pressures (to 1000 atmos- 
pheres), at different temperatures (25 to 55°C) and high shear 
rates were those of Appeldoorn, et al (1) and Philippoff (11). 
They employed the principle of the oscillating quartz crystal 
viscometer introduced by Mason (10) . More recently, Rein 
(13, 14) extended the techniques used by Philippoff to develop 
a viscometer for pressures to 8000 atmospheres. Rein subse- 
quently redesigned the entire apparatus to permit viscosity 
measurements to pressures approaching 15,000 atmospheres and 

temperatures to 250 “C and reduced shear rates on the order of 
8 

10 sec . In addition, the effect of pressure and temperature 
on the density of the fluids can also be measured. 

Experimental Equipment 

A flow diagram for the experimental equipment is 
shown in Figure 1. 

The high pressure cell, of shrink-fit duplex construc- 
tion with an 18 percent nickel maraging steel inner liner and 
an HY-140 outer shell, was designed conservatively for use at 
200,000 psi and 450 °F. The quartz crystal was mounted on the 
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Figure 1. Flow Diagram of High Pressure Equipment and Accessories, 
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inside of the top closure and was connected to two electrical 
leads which entered through steel-pipestone conical-shaped plugs 
in the top closure. 

A cell containing a coil of manganin wire (the change 
in resistance of the wire is a well-known function of pressure) 
was used for pressure measurements. The resistance of the manga- 
nin wire was measured with a Carey-Foster bridge. It is estimated 
that this combination measures pressures with an accuracy of 
better than ± 1/2 percent for pressures greater than 10,000 psi. 

A sheath thermocouple entered the cell through the 
bottom closure. It was intended to use this thermocouple to 
measure the fluid temperature inside the high pressure cell. 
However, this thermocouple failed to operate midway through the 
test program; subsequently temperatures measured with a thermo- 
couple on the cell wall consistently agreed with the temperature 
of the fluid in the cell to within + 1/2 ®F after a 15 minute 
time interval , based on previous observations . 

The temperature of the high pressure cell was maintained 
by immersing the cell in an oil bath which was heated by an 
electric heater. Temperatures of the fluid inside the cell were 
maintained to ± one percent. 

An intensifier was used to pressurize the cell. The 
intensifier could be recycled by appropriate use of valves 1 
and 2, shown in Figure 1, in order to attain the desired pressure 
in the cell. The principal disadvantage of this repetitive re- 
cycling is that it accelerates the wear in the packing and the 
fatigue life of the metal components. 

The intensifier was constructed so that the advancing 
piston forces oil from the low pressure side, at the juncture 
between the low- and high-pressure pistons, into a sight glass 
as shown in Figure 2. Measurement of the amount of oil displaced 
was converted into piston travel for calculation of the densities 
of the test fluid in the high pressure side. The sight gauge was 
calibrated in 1/8 inch increments which provided sufficient 
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sensitivity to allow density measurements to be made within 
± one percent, as confirmed by calibrations with di- (2-ethylhexyl) 
sebacate for which density determinations had been reported 
previously (2) . 

A liquid separator containing a piston and 0-rings pre- 
vented communication of the fluid being tested with the hydraulic 
oil in the pumping section. 

Valves 1 and 2, identified in Figure 1, consist of a 
section of U-shaped tubing which can be immersed in liquid 
nitrogen to freeze the test fluid; for this reason they are re- 
ferred to as frozen valves. They were introduced by Babb (3), 
and they have proved to be absolutely leak-proof. 

Figure 3 shows the electronic components used to measure 
the resistance of the oscillating quartz crystal as a function of 
frequency for each viscosity data point. 

A prime requirement for the oscillator is high stability. 
The manufacturer specifies less than 0.003 percent drift in 10 
minutes after warmup for this particular oscillator, which is 
adequate since the frequency will not change significantly during 
any measurement. 

This particular impedance bridge has the advantage that 
the undesirable' shunt effect of long cables is eliminated by the 
bridge design. The manufacturer claims an accuracy of at least 
+ one percent for resistance measurements with the bridge, and 
this accuracy can be generally attained. However, an accuracy 
of + 5 percent at the high end of the dial seems more reasonable 
because of difficulty in reading it precisely at the high end. 

Precision is the most important requirement for the 
counter used to determine the frequency of crystal oscillation 
because the calculated viscosities depend on a difference between 
frequency measurements. For fluids with low viscosities, such 
as lubricants at high temperatures, the resonant frequency in 
the test fluid differs from f^, the resonant frequency in an in- 
viscid fluid, by a small amount, 10-20 Hz. Therefore, if a 
counting or gate time of 10 seconds is used, the resulting 
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frequency measurements will be precise to ± 0.1 Hz. This pre- 
cision is necessary because viscosity calculations depend on the 
difference between the resonant frequency in the test fluid and 
the inviscid fluid, or f - f^. For more viscous fluids the 
difference between f and f is greater; in this case a gate time 
of one second was used and the resulting frequency measurements 
were precise to within + one Hz . 

The detector has a narrow band filter which is advanta- 
geous because it eliminated the effects of stray or undesired 
signals from the bridge balance. Incorporated in the detector 
is a local oscillator, which is driven at the same frequency as 
the impedance bridge. Since the local oscillator was used to 
drive the counter, correct frequency measurements were assured. 

The quartz crystal is the sensor for viscosity measure- 
ments. Preparation of the crystal was the same as described by 
Rein (14) except that the ring of silver paint at the crystal 
midpoint was omitted and the electrodes were separated by tape 
during silver deposition. When the tape was removed following 
deposition, the electrodes remained effectively separated. This 
technique of electrode separation is superior to the earlier 
technique because it eliminated unnecessary crystal abrasion. 

The qudrtz crystal must be securely mounted to insure 
reproducible results since the slightest slippage of the crystal 
can cause a significant change in the resonant frequency (as deter- 
mined by calibration with an inviscid fluid) which in turn will 
introduce an appreciable error in the calculated viscosities. 

Soldered electrical leads have been found to be unreliable 
mountings (15) because the leads often break off from the crystal. 
For this reason subsequent investigators have employed a spring- 
clamp mounting (1, 14). However, for this study a superior 
mounting was devised; it is shown in Figures 4 and 5. Its main 
advantage is that the crystal is supported at the crystal ends 
so that the electrical leads are not required to support the 
weight of the crystal. Further, reliable electrical contact 
to the crystal is attained by contacts in set screws which can 







11 



Figure 5. 


The 60 kHz Crystal in the Improved Mounting. 


! 

I 


12 

be advanced to maintain contact with the crystal as shown in 
Figure 4. 

Viscosity Measurements 

Prior to making viscosity measurements on the test fluids, 
a rather time-consuming calibration procedure was required. To 
facilitate this exposition and to demonstrate the methods for 
reducing the raw data, a brief review of the pertinent equations 
follows. 

Pertinent equations ; Viscosity measurements by the 
oscillating quartz crystal technique are interpreted in terms of 
a complex viscosity fh* = - ih 2 ) or complex shear modulus 

(G* = + iG 2 = iwn*) which is related to the mechanical imped- 
ance, + iXj^, of the oscillating rod by 

= Gj/ai = 2Xj,Rj,/(up (1) 

"2 = V“ ' (2) 

where P is the fluid density and u = 2irf is the angular frequency 

I 

of the rod. n-j^ (or G 2 ) relates to viscous behavior, while G^ 

(or 112 ) relates to elastic behavior of the fluid. 

and are related to the measured resonant frequency, 
f, in a given fluid and to the measured resistance at resonance, 

Rg, at the particular temperature and pressure, by 




(3) 

and 

' <*E - 

(4) 


where the subscript o refers to values obtained in an inviscid 
fluid. K and are related to crystal properties and can be 
calculated (1, 15) or obtained from calibrations in Newtonian 
fluids of known viscosity and density (1, 4, 10, 11, 14, 15). 

Manipulation of Equations 1, 3 and 4 shows that 
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Rg = (a)pn/2)^/^ (5) 

f = - Kj (topn/2)^/^ (6) 

for Newtonian fluids (for which = Rj^) • Thus, it is apparent 

that is the slope, and R^^ the intercept, of the R^, versus 

(o)pn/2)^^^ curve. Likewise, is minus the slope, and f the 

^ 1/2 ® 

intercept, of the f versus ((Dpri/2) ' curve. Therefore, for 

different Newtonian fluids with known viscosity and density, plots 
1/2 

of (wpri/2) ' versus resonant frequency and resistance give the 
desired calibration constants. 

In addition, and are related to physical properties 
of the crystal by the following equations (1, 15). 

Kf = (TTp^)"^(r"^ + h"^) (7) 


K = 4ttL K, 
r c f 


L = (2irAfG ) 
c c 


where r = the crystal radius 

h = the length of the crystal 
P^ = the crystal density 

= the conductance at resonance or the reciprocal of 
the resistance at resonance 

Af = the difference between the two frequencies at which 
the conductance is 1/2 the minimum value. 

Calibrations based on Equations 5 and 6 are preferred 
to use of Equations 7, 8 and 9 because Equations 5 and 6 are more 
directly related to the actual measurements. K^, and 
have been reported to be independent of pressure (1, 14). Also, 
a change in f^ with pressure has been reported (13) . Therefore, 
calibrations were obtained at all temperatures and pressures at 
which lubricant viscosities were to be measured in order to 
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confirm the effect of pressure on f^, to detemine the effect of 

temperature on f , and to validate the invariance of K , K^. and 

o r r 

with pressure. Since the oscillating quartz crystal visco- 
meter technique had not been used previously under the combined 
conditions of high t^perature and high pressure, these calibra- 
tions served two additional purposes : to determine whether the 

simultaneous application of high temperature and high pressure 
had any effect on the pertinent parameters or on the stability 
of the crystal mounting. Although no effect on the pertinent 
parameters, per se, was detectable, some erratic behavior was 
encountered in attempts to make measurements at the highest 
temperature (482°F) which might be attributable to the delicacy 
of the crystal mounting. 

Values of the viscosity and elastic modulus obtained 
at discrete frequencies and at various pressures and temperatures 
can be converted to equivalent values at a single temperature and 
pressure and over a range of frequencies by use of the concept of 
reduced variables (1, 8, 11). As applied to this investigation, 
the reduced variables concept states that increasing pressure, or 
decreasing temperature, has the same effect on viscosity or elastic 
modulus as increasing shear rate. 

The reduced variables are defined by (1) : 

<«> 

®TP = <V’'><Po/'’><V''o>p'V''o>T <“> 

=lr = 

''ir “ ni(nynT)p(n„/np>T (13) 

where the subscripts T and P refer to the temperature and pressure 
of measurement and the subscript o refers to values at the common 
temperature and pressure to which all measurements are being 
reduced. The viscosities r)^ and rip are those measured in the 

limit of zero shear rate. 
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Calibrations ; The fluids selected for calibration pur- 
poses were methylcyclohexane , n-pentane and di- (2-ethylhexyl) 
sebacate since their viscosities and densities over the range of 
pressures and temperature of interest were available. 

The first step in the calibration procedure was to 
measure the resistance of the quartz crystal as a function of 
the frequency of the crystal when immersed in an inviscid fluid 
(either a vacuum or air at ambient conditions) . Similar measure- 
ments were then made with the crystal immersed in a given fluid 
at a particular temperature and pressure. These measurements 
were repeated for each fluid at the other prescribed levels of 
temperature and pressure pertinent to this study. The data on 
resistance versus frequency for each fluid at each level of 
temperature and pressure were subjected to a least squares, 
quadratic (or cubic, where necessary) fit to derive an equation 
which was then differentiated to obtain the minimum point. The 
frequency at which the resistance was found to be a minimum was 
taken to be the resonant frequency at that particular temperature 
and pressure for the given fluid. 

It should be understood that this same procedure had to 
be exercised for each test fluid as well as each calibration 
fluid. A typical plot is shown in Figure 6 for one of the test 
fluids, FN-3158, at 210°F and 20,200 psi, which is quite repre- 
sentative of the data obtained at all levels of temperature and 
pressure for both the calibration and test fluids. It should be 
noted that the abscissa in Figure 6 is given in terms of 
(f - 19,900) rather than f. The only purpose in using this shift 
was to spread the data along the horizontal axis to facilitate 
visual examination. 

After measurements on a fluid were completed, the cell 
and high pressure lines were flushed several times with petroleum 
ether to remove any remaining test fluid. The system was then 
put under vacuum for several hours. Next, the resonant frequency 
and resistance were checked in vacuum to determine if the cleaning 
had been effective. If the resonant frequency or resistance 
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SHIFTED FREQUENCY = (f - 19900) Hz 

Figure 6. Typical Data for Determining Resonant Frequency. (Data for FN-3158 
at 210°F and 20,200 psi.) 
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deviated significantly from the original values in a vacuum (or 
air) , the cell was again flushed with petroleum ether and the 
test repeated. Measurements on the next fluid began only after 
it was obvious that the system was clean. 

Having obtained the relationship between resonant fre- 
quency, f, and resistance, R^, for the calibration fluids for 
which the viscosities and densities at various pressures and 
temperatures were known, it was possible to derive the calibra- 
tion parameters, K^, K^, f^ and as defined by Equations 1 

through 4. It is evident from Equation 5 that a plot of R_ 

1/2 ^ 
versus (topri/2) ' should yield a straight line having a slope 

of K and an intercept of R„ for each pressure and temperature 

^ 1/2 
level. Similarly, by Equation 6 , a plot of f versus (a)pri/2) ' 

should yield a straight line having a slope of -K^ and an inter- 
cept of f^ for each pressure and temperature level. The result- 
ing values for these calibration parameters, K , K-, f . and R_ 

3T r o Jbo 

then reveal to what extent they are dependent on pressure and/or 
temperature. 

The foregoing plotting technique was used for all cali- 
bration runs with the 20 kHz crystal. However, this method 
could be used to calibrate the 60 kHz crystal only at atmospheric 
pressure. At higher pressures, the inherent surface roughness 
of the 60 kHz crystal introduces intolerable uncertainties in 
the prediction of the calibration parameters (4) . Measurements 

obtained with calibration fluids such as n-pentane or methyl- 

1/2 

cyclohexane fall in the lower range of values for (o)pri/2) ' 
where surface roughness causes adverse effects. Since other 
fluids having known density and viscosity variations with pres- 
sure over the required range of high pressures are not available 
as calibration standards, the less preferred technique of using 
Equations 7 through 9 was employed to compute the calibration 
parameters at elevated pressures for the 60 kHz crystal. 

Although methylcyclohexane and n-pentane were suitable 
for all calibration runs on the 20 kHz crystal, di- (2~ethylhexyl) 
sebacate could be used only over a limited range of pressures 
and temperatures as a calibration medium. As will be discussed 
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later, it demonstrates viscoelastic behavior at higher pressures 
and lower temperatures (below 300°F) . 

The required data on viscosity and density as functions 
of temperature and pressiare for di- (2-ethylhexyl) sebacate were 
obtained from the ASME studies (2). Viscosity versus pressure 
data at 86 and 167®P from Bridgman (5) were interpolated to obtain 
viscosity-pressure information at 100 °F for n-pentane and methyl- 
cyclohexane. These same data were also extrapolated to 210°F 
and 300 “F to obtain the required viscosity-pressure information. 
Density-pressure information for n-pentane at 100 and 210 “F was 
obtained by interpolation and extrapolation of data reported by 
Bridgman (6) . 

The remaining density information for use in calibra- 
tions had to be calculated because experimental data are not 
available. Density calculations were based on a corresponding 
states scheme (12) . 

Figure 7 is representative of calibrations to obtain 

K and R„ for the 20 kHz viscometer. Since it includes data 
r Eo 

for all temperature measurements, the good correlations suggest 
the absence of a temperature effect of K^, which confirms con- 
clusions reached by others (1). was found to be independent 

of pressure and temperature in accordance with other investiga- 
tions (1). However, in contrast to the results reported by 
others (1, 14), K was found to depend on pressure, as indi- 
cated in Figure 7. A close inspection of the data given by 
Appledoorn, et al . (1) indicates a small change in with 
pressure. Since the magnitude of this variation would have a 
neglible effect on their results, it might have been overlooked. 
Uncertainties in data may have masked any changes in with 
pressure in Rein's earlier work (14). 

As explained before, the calibration parameters for 
the 60 kHz crystal at elevated pressures had to be calculated 
from Equations 7 through 9. The calculated values for for 
the 60 kHz crystal likewise showed a variation with pressure. 
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The dependency of on pressure for both crystals is summarized 
in Figure 8 . 

Figure 9 is representative of calibrations to obtain 
and f^ for the 20 kHz crystal. As indicated by the random 
variation in with pressure, it was concluded that was in- 
dependent of pressure which is consistent with previous observa- 
tions (1, 14). was also found to be independent of temperature, 

again in agreement with conclusions by others (1) . The change in 
f^ with pressure at 100 ®F as found in this study is in accord — 
at least within experimental error — with others (13) . Since f^ 
also varies with temperature, it had to be determined at every 
temperature and pressure for which viscosity measurements were 
to be made on the four test fluids. 

Table I summarizes the principal outputs from the 
calibration runs. 

Density Determinations 

It is apparent from Equations 1 and 2 that the density 
of the fluid must be known in order to compute the viscosities. 
Since density data at the different levels of temperature and 
pressure were not available for the four test fluids, they had 
to be measured in the course of this contract. The procedure 
consisted of introducing the test fluid into the pressure cell 
by means of the intensifier (refer to Figures 1 and 2) . The 
cell is maintained at the desired temperature level. The mass 
of test fluid introduced into the cell to attain a prescribed 
pressure was determined from the intensifier piston travel which, 
in turn, was indicated by the amount of pressurizing oil dis- 
placed into the sight glass. The experimental procedure and 
data reduction are elucidated in the following derivation of 
the pertinent equations . 

The mass of fluid introduced into the pressure cell, 
dmceii f is equal to the mass of fluid displaced by the intensi- 
fier piston, -dm, . Thus, 

T* on ' 


dm 


‘cell 


= -dm . . 

xnten 


= P 


inten 


Adil 


(14) 
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Figure 9. Typical Frequency Calibration Data. (Data from Run 20 kHz-: 
300°F.) 
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TABLE I 


CALIBRATION PARAMETERS FOR EACH DATA RUN 


Crystal* 

Parameter 

Value 

Comments 

20 kHz-11 

^f 

0.32 

Independent of temperature 
and pressure 


K 

r 


Independent of temperature, 
dependent on pressure as 
shown in Figure 8 


f 

o 


Empirically determined for 
each temperature and pres- 
sure of measurement 


^Eo 

31.6 kJ^ 

Independent of temperature 
and pressure 

20 kHz-111 

^f 

0.30 

Independent of temperature 
and pressure 


K 

r 


Independent of temperature^ 
dependent on pressure as 
shown in Figure 8 


^o 


Empirically determined for 
each temperature and pres- 
sure of measurement 


^Eo 

30.4 kn 

Independent of temperature 
and pressure 

60kHz** 

^f 

0.32 

Independent of pressure 


K 

r 


Dependent on pressure as 
shown in Figure 8 


^o 


Determined from correlation 
in Reference 13 


^Eo 

Ofi 

Independent of pressure 


*20 kHz-11 and -III designate re-mounting of the crystal. 
Calibrations for runs -0 and -I were not obtained because of high 
pressure equipment failures during these runs. 

**The temperature was constant during all of these runs so 
no variation in parameters with temperature could be observed. 
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where A is the cross sectional area of the high pressure piston, 
Pinten density of the fluid in the intensifier at the 

existing pressure and temperature of the intensifier, and d5, is 
the piston movement during pressurization. 

The total mass in the cell at any pressure is the mass 
that was present at atmospheric pressure, 
added during pressurization. Thus, 


m^, and the amount 


m T , = ra + 
cell o 




dm = m + 
cell o 




p . . Ad)l 
^inten 


(15) 


Division of Equation 15 by the volume of the cell, V, a constant, 
gives 


P = Pcell ' Po j Pinten 

= Po + (A/V) j dP 


(16) 


where p = the fluid density at a given temperature and elevated 
pressure 

pQ = the fluid density at the given temperature and atmos- 
pheric pressure 

A = the cross sectional area of the high pressure inten- 
sifier piston 

V = the volume of the high pressure cell 
(3)l/3P),j, = the slope of the piston displacement-pressure curve. 

The restriction of constant temperature was satisfied 
by carrying out the compression slowly. 

The cell volume, V, in Equation 16 can be approximated 
closely by 


V = Total internal volume of cell - Crystal volume - 
Volume of viscometer support 

From the measured dimensions of these quantities, the cell volume, 

3 

V, was calculated to be 60.9 in . 
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The high pressure piston had a diameter of 0.625 inches 
■ 1 ::; 2 
(thu|, a cross sectional area of 0.307 in ) during all measure- 
ments, except for the runs with Humble FN-3158 fluid. In this 

latter case an oversize piston having a diameter of 0.660 inches 

2 

(cross sectional area of 0.342 in ) was used. 

In addition to the values of V and A obtained above, 
it is necessary to know a function of pressure at the 

intensifier temperature. Based on the material balance of Equa- 
tion 14, and restricting it to the case of equal cell and intensi- 
fier temperatures, the following finite difference equation was 
used to determine ^ function of pressure; 


Pn+i - Pn = [1 + (A/2V)AA^]/[1 


(A/2V) AA^] 


(17) 


M represents the piston displacement during a given pressure 
change and the subscript T indicates that the intensifier and 
cell must be at the same temperature. Densities of the test fluids 
were measured at temperatures above room temperature and the in- 
tensifier was always at room temperature. Thus no measurements 
were made in which the cell and intensifier were actually at the 
same temperature. Therefore, AJl^, at given constant pressures, 
was obtained by extrapolation of the measured values of AS, at 
the various temperatures to the intensifier temperature. 

The foregoing technique for measuring densities was 
validated by measuring the densities of one of the calibration 
fluids, di- (2-ethylhexyl) sebacate, for which density data were 
available in the literature. Table II compares these values 
and demonstrates an agreement within + one percent. 
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TABLE II 

A COMPARISON OF DENSITY VALUES FOR 
DI- (2-ETHYLHEXYL) SEBACATE 


Temp. 

(°F) 

Press . 
(1000 psi) 

p (This Project) 
(gm/cc) 

p (Reference 2) 
(gm/cc) 

100 

10 

0.939 

0.941 


20 

0.968 

0.965 


40 

1.008 

1.007 


60 

1.034 

1.037 


80 

1.057 

1.064 

210 

10 

0.903 

0.906 


20 

0.938 

0.934 


40 

0.985 

0.982 


60 

1.014 

1.013 


80 

1.037 

1.064 


100 

1.067 

1.064 

300* 

10 

0.870 

0.879 


20 

C).908 

0.913 


40 

0.960 

0.963 


60 

0.989 

0.997 

482** 

10 

0.813 

0.832 


20 

0.866 

0.874 


40 

0.935 

0.928 


60 

0.968 

0.965 


80 

0.986 

0.996 


100 

1.013 

1.023 


115 

1.054 

1.039 


*Density values for Reference 2 were obtained by inter- 
polation of known p (T) data at each pressure. 

**Density values for Reference 2 were obtained by extra- 
polation of known p (T) to 482°F at each pressure. 






RESULTS 


Having determined the calibration parameters, f^, K^, 
and RgQf for the viscometer, measurements of frequency and 
resistance of the crystal in the four test fluids as well as the 
density of the fluids are required at the various levels of tem- 
perature and pressure in order to calculate the viscosities from 
Equations 1 and 2. The experimental procedures are identical to 
the calibration exercises described previously. This section 
presents the results of the density and viscosity determinations. 

In addition (although not part of the contractural re- 
quirements) the observations of viscoelastic behavior are discussed 
briefly as relevant information. Indications of viscoelastic 
behavior in Mobil XRM-109 and Humble FN-3158 plus 10 volume per- 
cent Kendall 0839 were found at 100°F, but Sufficient data were 
not obtained to be conclusive. However, to illustrate the analysis 
for viscoelasticity, the results obtained with one of the calibra- 
tion fluids, di- (2-ethylhexyl) sebacate, are presented. 

Density 

The densities of the four test fluids as functions of 
pressure at different temperature levels are siammarized in Figures 
10, 11 and 12. The densities at atmospheric pressure for various 
temperatures were obtained by interpolation and extrapolation of 
data supplied by NASA (9) as shown in Figure 13. Tabulated density 
data are given in Appendix A. 

As could be expected, the density-pressure curves for 
Humble FN-3158 and the mixture FN-3158 plus 10 volume percent 
Kendall 0839 were found to be essentially identical. The density 
of each of these components is similar at atmospheric pressure, 
and apparently the presence of 10 volume percent of Kendall 0839 
does not introduce any deviations from additive volume behavior. 
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DENSITY, 6M/CM 



PRESSURE, 1000 PSI 

Figure 10. The Effects of Pressure and Tempera 
ture on the Densities of Humble 
FN-3158 and Humble FN-3158 + 10 vol 
% Kendal 0839. 



20 30 

PRESSURE, 1000 PSI 


40 


fects of Pressure and Temperature on 
nsity of Mobil XRM-109. 



DENSITY, GM/CM 
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Figure 12. The Effects of Pressure and Temperature on the 
Density of DuPont Krytox 143-AB. 




DENSITY, GM/CM 
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10^ X T/T, 

Figure 13. Density of Test Fluids at Atmospheric 
Pressure. (Data from Reference 9.) 


DENSITY, GM/CM 




32 


Viscosity 

Figures 14, 15, 16 and 17 summarize the results of the 
viscosity measurements for the test fluids as functions of 
pressure at different temperatures. Tabulated data are given in 
Appendix B. The viscosity data are estimated to be accurate 
with ± 15 percent. It will be noted in Figure 17 that viscosity 
data for DuPont Krytox 143-AB are given only at 100°F. Although 
data were taken at higher temperatures, they were not considered 
sufficiently reliable to report. Upon completion of the higher 
temperature runs, the pressure cell was opened. Examination of 
the crystal mounting revealed that it was not true. Time was not 
available to repeat these runs at the higher temperatures, par- 
ticularly since the GR 1310-A oscillator was in need of maintenance 
which required a down-time of several weeks. 

For the sake of comparison, the viscosities of the four 
test fluids at atmospheric pressure at temperatures of 74.3°, 

112°, 130° and 190°F were determined by means of a Canon-Fenske 
capillary viscometer. These data are summarized in Figure 18. 

No attempt was made to extrapolate these data beyond 210°F because 
of the uncertainties involved. The capillary values are shown on 
the zero gauge pressure ordinate for the 100° and 210°F isotherms 
in Figure 14 through 17 . The discrepancies between the two methods 
of measurement are within the estimated accuracy of ± 15 percent 
for the oscillating quartz crystal viscometer. 

Viscoelastic Behavior 

Although the work statement for this contract did not 
include determination of the viscoelastic effect, the fact that 
two of the test fluids, Mobil XRM-109 and Humble FN-3158 plus 
10 volume percent Kendall 0839, gave indications of departure 
from Newtonian behavior merits mention. 

The components of the complex viscosity of a fluid 
are given by 


n* = - iri2 


(18) 



VISCOSITY 












Figure 17. The Effect of Pressure on the Viscosity of 
Dupont Krytox 143-AB at 100*F. 




VISCOSITY, POISE 
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RECIPROCAL TEMPERATURE, 10^ x 1/T, °R“^ 


Figure 18. Viscosity as Determined from Capillary 
Viscometer at Atmospheric Pressure, 
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and the components of the complex shear modulus by 

G* = + iG2 = iwn* (19) 

Both are related to the mechanical reactance, Xj^, and resistance, 
of the oscillating crystal by 


= G2/0) = 

2XmV“‘= 

(1) 

= G^/o) = 


(2) 


where P is the fluid density and w = 2irf is the angular frequency 
of the crystal, which is assumed to be equal to the shear rate. 
The mechanical impedance of the crystal, Zj^, is given by + iX^^ 
where 


II 

- Rtp ) /K 
Eo ' r 

(4) 

o 

II 

- f)/Kf 

(3) 

For a Newtonian fluid, ri* = n^. 

since. r\^ = 

= 0 in Equation 18 , 


which requires that Rj^ = X^^ according to Equation 2 . However , 

if ri» is non-zero, then R,, 4 X„ and viscoelastic behavior is 
2 M M 

indicated. 

A characteristic of the oscillating quartz viscometer 
is that the difference (R^ - in Equation 4 is much larger 

than the difference (f^ - f) in Equation 3. For this reason, 
errors in measurements of resistance are less serious than errors 
in measurements of frequency for computing values of Rj^ and Xj^, 
respectively. As the temperature of measurement increases and 
simultaneously the viscosity decreases, f approaches f^ much 
faster than R^ approaches Consequently, a loss in signi- 

ficant figures for the difference (f^ - f) results. Because of 
this . uncertainty , whenever the calculated values of Rj^ differed 
significantly from X^^ for measurements at a particular temperature 
and pressure, whereas measurements on this fluid consistently 
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resulted in Rj^ being equal to at other pressures and tempera- 
tures (for which a difference in and would have been more 
plausible) the calculation of viscosity, was based entirely 
on the more reliable resistance measurements. In this case 
Equation 1 becomes 

= 62 /“ = 2 Rj^V“P (for ~ ^ 20 ) 

The viscosities which were calculated solely from Equation 20 
are indicated by enclosing the data points in parentheses for 
the 300®F isotherm for Mobil XRM-1 09 in Figure 15 and the 210®F 
and 300®F isotherms for Humble FN-3158. Aside from these excep- 
tions, both frequency and resistance measurements were used to 
compute the viscosities by means of Equations 1 and 2. 

For all fluids at 100®F the difference (f^ - f) retained 
sufficient significant figures so that reliance could be placed 
on the frequency, as well as the resistance, measurements. In 
other words, at 100“F equal confidence could be given to the 
values of and ^2 as computed from Equations 1 and 2. The 
components of complex viscosity, complex shear 

modulus, + iG 2 , are summarized in Table III for Humble FN-3158 
plus 10 volume percent Kendall 0839 and for Mobil XRM-109. For 
the other two fluids. Humble FN-3158 and DuPont Krytox 143-AB, 

02 was found to be zero. The shear rates in the fluid were 
obtained by making the usual assumption that they are equal to 
the angular frequency of the crystal, w, which is defined as 

OJ = 2iTf. 

By resorting to the use of the concept of reduced 
variables, as defined by Equations 10 through 13, the onset of 
viscoelastic behavior is predicted at 100 °F and atmospheric pres- 
sure for the following reduced angular frequencies (or shear 
rates) : 



VISCOELASTIC BEHAVIOR OF TEST FLUIDS AT 100 °F 
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Reduced 

Fluid Angular Frequencies 

g 

Humble FN-3158 plus 4.7 x 10 radians/sec 

10 volume percent 

Kendall 0839 

Mobil XRM-109 5.9 X 10^ radians/sec 

The frequencies listed above are at best approximate because 
viscosities at zero shear rate as a function of pressure and 
temperature, which are required to compute the reduced variables, 
were not available. Also, until more experimental data are ob- 
tained over a wider range of conditions, the conclusion of visco- 
elastic behavior for these two fluids at 100°F should be regarded 
as tentative. 

The difficulty with the present design of the high pres- 
sure equipment is that it was not possible to make measurements 
above 40,000 psi because of the excessive viscosity of these 
fluids. Furthermore, at higher pressures, a deterioration of 
sensitivity in measurement was encountered. The principal 
advantage of going to higher pressures is that it would permit 
simulation of higher shear rates at atmospheric pressure by 
means of the concept of reduced variables. Further alleviation 
of the present' limitations may be forthcoming with the 100 kHz 
crystal, which is now being prepared for preliminary testing 
in the viscometer. An alternative may be the use of another 
quartz crystal technique reported by Barlow and Lamb (8) , for 
which shear rates of at least two orders of magnitude higher 
than the present oscillating quartz crystal can be simulated. 

During the calibrations of the viscometer with di- (2- 
ethylhexyl) sebacate, it was observed that the measured fre- 
quency was significantly higher than would be expected for a 

1/2 

Newtonian fluid for given values of ({opri/2) ' at the higher 
pressures. The anomalous behavior was consistently observed 
for both crystals and for the several remountings of the 20 kHz 
crystal that were required during the course of the experimental 
work; these observations confirmed the existence of viscoelastic 
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behavior. The temperatures and pressures at which these devia- 
tions from Newtonian behavior were observed are summarized in 
Table IV. 


TABLE IV 

CONDITIONS FOR VISCOELASTIC BEHAVIOR OP 
DI-(2-ETHYLHEXYL) SEBACATE 


Crystal 

Temperature 

Pressure 


(®F) 

(1000 psi) 

60 kHz 

100 

10 and 18 

20 kHz 

100 

10, 20, 40 and 60 

20 kHz 

210 

20, 40 and 60 


The concept of reduced variables was again employed 
to convert these observations with the two crystals at various 
temperatures and pressures to a common basis of 100°F and atmos- 
pheric pressure. Figure 19 shows the reduced elastic (storage) 
modulus, as a function of reduced frequency. The results 

of previously reported observations on the viscoelastic behavior 
of di- (2-ethylhexyl) sebacate at much higher frequencies (8) 
are also shown. A least squares fit of the data from both of 
these investigations is given by the equation 

Gir = 0.189 (2i) 

with a correlation coefficient of 0.992. (Note that the data 
points taken from Reference 8 and showh on Figure 19 were read 
from a smoothed curve of R^^ and Xj^ versus reduced frequency; 
thus they exhibit less scatter than the raw data from the pres- 
ent investigation.) Even though there is no overlap of the data 
between the two investigations, the high correlation coefficient 
indicates that both sets of results are in basic agreement. 




Figure 19. The Storage Modulus of di- (2-ethylhexyl) Sebacate 
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Furthermore, this agreement serves as a confirmatipn of the con- 
cept of reduced variables. 

The correlation shown in Figure 19 was used to calculate 
a shear relaxation spectrum for di- (2-ethylhexyl) sebacate by 
means of the defining Equation (7) ; 

H (at T = 1/w) = 2 /tt Im (at (22) 

where H is the relaxation spectrum 

Im signifies the imaginary part of evaluated at 
This spectrum is shown in Figure 20; a relatively wide distribu- 
tion of relaxation times is indicated. 






CONCLUSIONS 


The principal conclusions resulting from this investi- 
gation are; 

1. The present oscillating quartz crystal viscometer can be used 
to measure viscosities at elevated temperatures and pressures 
and discrete shear rates to an accuracy of at least ± 15 
percent. 

2. The concept of reduced variables in combination with viscosity 
measurements at high pressures to simulate high shear rates 

at atmospheric pressure has been reconfirmed. 

3. Viscoelastic behavior is manifest at atmospheric pressure and 
100®F for Mobil XRM-109 fluid at an angular frequency of about 
5.9 X 10^ radians per second and for Humble FN-3158 plus 10 
volume percent Kendall 0839 at an angular frequency of about 
4.7 X 10^ radians per second. 

4. Previous observations that the calibration parameters for the 

crystal viscometer, and are independent of both pres- 

sure and tejnperature are confirmed. The dependency of the 
calibration parameter, f^, on both pressure and temperature 

is also consistent with earlier studies. However, the observed 
change in the calibration parameter, K^, with pressure has 
not been reported previously. 

5. The viscoelastic behavior of one of the calibration fluids, 
di- (2-ethylhexyl) sebacate, as measured by the reduced elastic 
(storage) modulus, corresponds closely with data reported by 
Lamb in 1965 using a different technique. The calculated 
relaxation spectrum shows no maxima and is characteristic of 

a broad distribution in relaxation times. 

6. A thorough investigation of the possible viscoelastic effects 

in fluids would require three experimental techniques: (1) 

The effect of pressure and temperature on viscosity at zero 
(or maybe very low) shear rates such as the falling-ball 
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method; (2) the effect of pressure, temperature and variable 
shear rates up to several million Hertz by the oscillating 
quartz crystal technique; and (3) the effect on viscosity of 
very high shear rates, at least up to 10^^ Hertz by the 
quartz technique utilized by Barlow and Lamb. 



LIMITATIONS AND REMEDIES 


During this study several deficiencies in the design 
of the oscillating quartz crystal viscometer became evident which 
limited the range of feasible measurements. These shortcomings 
have since been corrected, or can be corrected, with a modest 
cimount of development work. The fact that this particular equip- 
ment was entirely new and had not been used previous to the ini- 
tiation of this study probably accounts for several of the experi- 
mental limitations uncovered during the earlier phases of the work. 
The major deficiency in the present design is the difficulty in 
handling very viscous fluids (10 poise) such as were encountered 
in this study. 

Sensitivity 

As the viscosity of the test fluid increases, the sensi- 
tivity of the measuring equipment decreases. The sensitivity, 
which is given by the ratio of the change in bridge output to the 
change in the resistance, is compromised as the resistance in- 
creases with in'creasing viscosity. For example, in this study 
it was very difficult to detect changes in resistances of less 
than 5 percent. Also, the dial on the bridge which indicates 
the resistance values has decreased resolution at the high end. 

The limitation on dial resolution has been circumvented 
by inserting a 1.500 Mfi resistor in parallel with the crystal, 
thus lowering the measured resistance and thereby keeping the 
measured resistance within a readable range on the dial. This 
technique was used throughout the present investigation. Of 
course, this adaptation further reduced the sensitivity of the 
bridge output to crystal resistance changes, but the deleterious 
effect was more-than-compensated-for by the increased ability to 
read the dial. Also, with the parallel resistor, the actual 
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crystal resistance had to be calculated from the resistance 
relation for parallel resistors. 

At present, further amplification of the bridge output 
signal, to increase measurement sensitivity effectively, is 
limited by noise associated with the signal. Decreasing the 
noise to permit greater amplification seems to be the most 
promising way to increase measurement sensitivity. Sources 
of noise, in order of decreasing probability of causing problems/ 
are: (1) pickup by the leads from the crystal to the measuring 

equipment, (2) the oscillator and (3) mismounting of the crystal. 

High Pressure Cell Design 

In the present design, the test fluid is pressurized in 
the intensifier and is transferred into the high pressure cell 
containing the viscometer through a transfer line having an in- 
ternal diameter of 0.025 inches. At higher pressures (above 
40,000 psi) the viscous test fluids behave virtually as solids 
in the transfer line and it is not possible to pressurize the 
cell. Some preliminary attempts were made to heat the intensi- 
fier and transfer lines, but this remedy introduced more problems 
(with respect to the operation of the intensifier and the man- 
ganin coil) than it solved. The most promising solution is to 
enclose the crystal viscometer in a bellows. The test fluid 
would then be charged directly into the bellows chamber at essen- 
tially atmospheric pressure. It would then be brought to the 
desired pressure by using the intensifier and a low-viscosity 
fluid such as pentane to transmit pressure across the bellows 
to the test fluid. A schematic of this proposed design is shown 
in Figure 21. Another advantage for the bellows chamber is that 
it would reduce greatly the amount of test fluid needed for 
measurements . 

Electrode Integrity 

The silver electrodes on the quartz crystal are vacuum 
deposited. In the earlier stages of this investigation, it was 
observed that the silver tended to flake-off after pressurization 
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ire 21. A Schematic Diagram of a Bellov 
for Containing the Test Fluid j 
Pressure Cell. 
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at temperatures cJ?ove 300®F, thus breaking the electrical circuit. 
This problem now appears to have been resolved by maintaining 
closer control of the thickness of the silver deposit. 



RECOMMENDATIONS 


1. The viscoelastic behavior of Mobil XRM-109 and the various 
mixtures of Hximble FN-3158 with Kendall 0839 should be 
investigated further at higher pressures by the quartz crystal 
viscometer technique, particularly with a 100 kHz crystal. 
Although most of the measurements should be made at 100 °F, 
for which viscoelastic behavior has been detected, measure- 
ments at higher temperatures would be desirable for conclu- 
sive confirmation of the data. 

2. Additional measurements should be made on di- (2-ethylhexyl) 
sebacate at elevated pressures with the crystal viscometer 

to provide data for "filling-the-gap" in the relaxation spec- 
trum between data from this investigation and a previous 
investigation. 

3. The combination of the crystal viscometer and high pressure 
should be used for measurements on any fluid for which 
viscoelastic behavior is expected. 
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APPENDIX A 


TABULATION OF CALCULATED DENSITIES 
AND PISTON DISPLACEMENTS 



Fluid Temperature 

Pressure 

Displacement 

Density 

(°F) 

(1000 psi) 

(in) 

(gm/cc) 

FN-3158 100* 

atm 

mmimm 

0.873 


10.0 

6.500 

0.905 


20.1 

10.937 

0.932 


39.8 

17.812 

0.970 

210* 

atm 

— — 

0.836 


10.0 

7.312 

0.872 


19.9 

12.750 

0.902 


40.4 

20.375 

0.945 

300* 

atm 


0.803 


10.0 

9.000 

0.846 


19.9 

14.812 

0.882 


39.5 

22.812 

0.925 

480 

atm 

— 

0.739 


10.0 

15.750 

0.817 


19.9 

23.625 

0.858 


39.7 

33.500 

0.895 

FN-3158+ 100 

atm 

— 

0.873 

10 vol % 0839 

9.8 

7.125 

0.905 


19.9 

12.375 

0.932 

* 

38.2 

19.500 

0.970 

210 

atm 

— 

0.836 


10.1 

8.125 

0.871 


20.0 

13.937 

0.902 


40.0 

22.312 

0.945 

300 

atm 

— 

0.803 


9.9 

9.875 

0.846 


20.1 

16.625 

0.882 


40.0 

25.250 


482 

atm 



0.739 


10.0 

17.375 

0.817 


20.1 

26.375 

0.858 


*Measurements at these temperatures were made using an 
oversize (0.660 dia.) piston. (Regular piston is 0.625 dia.) 
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Fluid 

Temperature 

Pressure 

Displacement 

Density 


(°F) 

(1000 psi) 

(in) 

(gm/cc) 


XRM-109 


Krytox 


100 

atm 

— 

0.839 


10.0 

7.250 

0.871 


20.0 

12.500 

0.899 


36.3 

18.825 

0.937 

210 

atm 

— 

0.806 


10.0 

8.125 

0.842 


20.0 

14.187 

0.872 

300 

atm 


0.778 


9.9 

9.750 

0.821 


19.8 

16.438 

0.855 

100 

atm 


1.88 


10.0 

10.000 

1.96 


19.9 

16.375 

2.04 


25.3 

19.250 

2.06 

210 

atm 

— 

1.76 


10.2 

12.750 

1.88 


19.7 

19.875 

1.97 


28.0 

24.125 

2.02 

300 

atm 



1.67 


10.1 

16.625 

1.82 


19.9 

24.750 

1.92 


23.8 

27.250 

1.95 


26.4 

28.500 

1.96 

437 

atm 


1.50 


10.0 

22.375 

1.73 


19.9 

32.000 

1.87 


21.9 

33.500 

1.89 





APPENDIX B 


TABULATED VISCOSITY DATA 



DATA FOR VISCOSITY CALCULATIONS AND 
CALCULATED VISCOSITIES 
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